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Abstract 

A new method of calculating deep inelastic structure functions for 
the nucleon in independent particle models is presented. The method is 
applied to the bag-like model and its predictions are compared with the 
parametrisation of Gluck, Reya and Vogt. 


1 Introduction 

Deep inelastic scattering provides worthy information about the structure of 
nucleons and nuclei. The papers of Gluck, Reya and Vogt ||, || have shown that 
quite good agreement with the data (and moreover some predictions) may be 
obtained assuming valence-like parton distributions at some small scale Q 2 = 
/P and evolving them using perturbative QCD. Nowadays it is not possible 
to calculate the initial distributions using directly QCD and we are left with 
various, less or more phenomenological, models. The idea of calculating leading 
twist parton distributions using quark models has over twenty years Q and 
there were many implementations of it using the bag model, soliton models, the 
colour-dielectric model and others ||. 

It is reasonable to assume that interactions between quarks may be approx¬ 
imately described by a mean field confining them. I fact, even one of the most 
prominent quark models, the bag model [[Tl] (in the static cavity approximation) 
may be treated in such a way. 

In the mean field approach quarks move independently and there arises the 
so called centre of mass problem, resulting from the simple observation, that 
independently moving particles cannot stand for a composite state of definite 
total momentum. The problem manifests itself in the fact, that quark distribu¬ 
tions, calculated in such a model, spread over all Bjorken x from —oo to +00 
instead of having support in the region from 0 to 1. There are many methods 
of curing this problem, but all of them have some drawbacks. For example 
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the method using Peierls-Yoccoz or Peierls-Thoules projection gives badly nor¬ 
malised quark distributions. In this paper we follow the idea of Szymacha Q 
that the independent particle model of a composite particle corresponds to the 
particle moving in an auxiliary external potential. 


2 Deeply inelastic nucleon structure functions 


Tensor W describing inclusive photon-nucleon scattering may be written as 

2p° 


47r 


J d 3 x j d 3 y^2(h\jrt(x)\n)(n\r{y)\h)e- iir( - s -y ) 2Tr6(p 0 + q 0 -P°), (1) 


where q° and q are the photon energy and momentum, |ro) are normalised to 
1 final states with energy P°, and |ft) is the initial (normalised to 1) nucleon 
state with energy p°. 

For | ft.) being an eigenstate of definite momentum p and spin s 


ft) 


1 

TPf 


h(p, s)), 


( 2 ) 


where (ft(p', s')\h(p, s)) = (27 T) 3 2p°5 a ' s 5 3 (p' — p) and (formally) V = f d 3 x = 
(27r) 3 ft 3 (0), 

we obtain usual form for IY A “' 

= <*k Mt (0)|n)<n| J"(0)|ft)(27r) 4 <5 4 (p + q - P n ). (3) 


If we approximate the real nucleon by independent particles (as in the bag 
model), then I ft) cannot be a state of definite momentum and (li) is much more 
useful then (|j). 

W may be expressed in terms of scalar functions Pi, F 2l gi and g 2 : 


= -{g^- 


q> x q v 


2 )F 1 + (p»-™q»)(p''-™f 


F 2 
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( 4 ) 


where m is the nucleon mass. For specific choices of p and v we obtain 

w tt =iv«= ( 5 ) 

= Pi, (6) 


where Q = yj-q 2 , x = Q 2 /2pq, n f = (p^ - ^q ^ /\jm 2 - and _L, U 

denote components perpendicular to n and q^. 
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3 “Naive” structure functions in the bag model 

As was already stated, independent particles cannot correspond to a state which 
is a momentum eigenstate. Following Szymacha et. al. j^] we will assume that 
they describe the nucleon with (yet unknown) wavefunction 4 'cm, bounded 
by an auxiliary external potential. Nevertheless our first step is to calculate 
structure functions of the bag model (in the static cavity approximation) without 
any improvements for the centre of mass motion. 

We assume that in DIS photon interacts only with a given quark from a 
nucleon, without any influence on other quarks, exiting it into the final state 
| n q ). Then Eq. ®) may be rewritten as (we add the subscript q to indicate that 
this W^ is calculated directly from the quark wave function) 

wr= TF f d 3 x f d 3 yYln q 

(g| (x)\n q )(n q \J v {y)\q)e~ lq ^ x ~' q) 2T:5{p 0 + q° - P°) = 

mJ26(p 0 + q°-P°) (7 ) 

n q 

J d 3 xe~^ q ~ Pn * )x $ q {x)Yu nq J d 3 yP < ' q - pn « )v u rlq Y'>Pq(y), 

where it was assumed that for higher excitation, final states |n 9 ) may be approx- 

_2^?^ X^ I 2 P X • • • • • 

imated by plane waves u nq e n i nq , and ip q is the initial wave function of 
the quark q. 

Since bispinor u Uq is uniquely described by corresponding momentum P nq 

and mass M = ~ (X 3 ) , we may replace E„, u n q u nq by Em I d 3 P 

- ^ P + M For large Q 2 we obtain 

(27t) 3 2VM 2 +P 2 

W r = T~Yl [ dae ~ wZa [ d 3 rip q (r}T^ilj q {f+e z a), (8) 

n m J J 

where p z = e q — mx( 1 + ^-) (e q is the initial quark energy), e z = and 

r “ = ^( 7 °- 7 z ), ( 9 ) 

= 7 °- 7 z - ( 10 ) 

We have to take into account that mass of the exited quark M is the same 
as the mass of quark in the ground state, hence there is only one term, with 
M = m q ~ 0, in the sum Em- 
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4 Taking into account the centre of mass motion 

The structure functions calculated in the previous section should be considered 
as structure functions for the nucleon in an artificial external potential. Photon 
interacting with such nucleon may excite it into higher eigenstate in the artificial 
potential, or excite its internal degrees of freedom, or do the both. Since the 
nucleon is a composite object, its coupling to the photon cannot be described 
by simple matrix 7 ^. We assume this coupling to be of the form 

f(q 2 ,P 2 ,Pq) 7 ", (11) 

where P M is momentum of the final state. 

In general this coupling may have more complicated spinor structure, but for 
large Q 2 other choices are effectively equivalent with ours provided we al¬ 
low f(q 2 ,P 2 ,Pq) to be dependent on the photon polarisation, i.e. we take 
f(q 2 ,P 2 ,Pq) = f A (q 2 ,P 2 ,Pq ) 1 where A = pu = tt, _LJ_. For large Q 2 , due 
to the energy conservation, we obtain Pq = —— y*—, so i n fact our formfactor 
depends only on the polarisation, Q 2 and M 2 = P 2 . For a given process (fixed 
q^) Q 2 remains unchanged and we may write W A “' for fiv = tt , _L_L in analogous 
to (||) form (hereafter A = fu>) 

= l/ A ( M2 )| 2 -T / da e~ ipta f d?r^> C m{?)T a ^ C m{t + e z a), ( 12 ) 

with p z = E — mx{ 1 + jjp-), and T u , r-*" 1 - given by Eqs. (||), (|To|). (fcM is the 
nucleon wave function, and E stands for its energy.) 

Introducing variable y = 1/(1 + 773 -), which by definition satisfies the con¬ 
dition 0 < y < 1 , the above equations may be rewritten in the form 

W A (x,Q 2 ) = f ^W A f {y,Q 2 )W A M (-,Q 2 ), (13) 

j 0 y y 

where 

W? ff (v) = (14) 

Wt f j{y) = \f^{y)\ 2 M{y), (15) 

(M(y) is a measure coming from replacing Ym by / M(y) = — 1 ) 

and 

W a m ^Q 2 ) = £ | dae-W-^J d^cM{^)T A M ^cM{?+S z a), (16) 
with (note difference with (§)) 

r cm = r hi* = 7° - i z - (it) 
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For a constant ^cm (corresponding to not moving, unbounded nucleon), we 
simply obtain 

W A M (^,Q 2 ) = ^2n6(E-m^). (18) 

In that case E = m and hence from ( [l3| ) 

W A = \w e A ff . (19) 

Comparing the above expression with (||), (;G[) we obtain 

Wt} f = — -2F U (20) 

JJ x 

Wt s i = 2 F u ( 21 ) 

This means that to obtain predictions for structure functions in our model, 
we have to find W A ff. Since for massless quarks W™ = 0 (the Callan—Gross 
relation), we are interested only in W^, which is directly related to the quark 
distribution function. To obtain we have to perform three steps: 

a) calculate the naive structure function W f j ~ x , using given quark wavefunction; 

b) calculate Wq^, using the centre of mass wavefunction (function of the 

nucleon in the auxiliary potential); 

c) solve the equation Eq. (fbl) for W A A using found W^m an d and taking 

w 3 - 1 - = W^. 

Of course as far as we don’t know the centre of mass wavefunction >1 'cm, 
this procedure cannot be performed. We will learn from the next section, how 
to find Wqm without calculating explicit form of ^cm- 

5 Calculating 

To obtain from Eq. (^G|), we have to calculate f d 3 r^cM{r)^'cM'^CM{r‘+ 
e z a). Since T^;^ = 7° — 7 Z , we may write 

J d 3 r^cM(r)T^ M ^CM{r + e z a) = 

f d 3 r^ C M{r)^CM{r + e z a) - Jd 3 r^ C M(r)^o: z 'i>cM(r + e z a). 

( 22 ) 

The term J d 3 f^CM{r+e z a) should correspond to overlapping function 
between translated states of the nucleon treated as a composite object. In our 
model the nucleon consist of three quarks with wavefunctions ip qi , so 

f d 3 f^ C M(f)^ cm{t + e z a) = (/ d 3 fip qi (r)^ qi (f + e z a)) ■ , , 

(f d 3 ft/j q 2 (r)^ q2 (r + e z a)) ■ (fd 3 fi/j q 3 (^i/j q3 (r + e z a)) . 
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To obtain analogous expression for / d 3 r^ ad’F CM{r + e z a) let us remind 
that for a Dirac particle, mean value of the matrix a may be interpreted as a 
mean velocity. It is thus natural to take 

&CM = -j^i e qi&qi + £ 92^*92 + £ 92 ^ 93 )) ( 24 ) 

where e qi is the energy of the quark qi , and E is the total energy of the system. 
For identical quark wavefunctions we obtain 

f d 3 r$>cM(r) 1 *u z '&CM{r + e z a) = 

(J d 3 fip qi (r}^a z ip qi (f+e z a)) (f d 3 rij} q 2 {f)^ qi (r + e») 2 . 

Equations (^3[) and ( |25| ) mean that we can calculate Wqm just using ij) q , without 
introducing '1 ; cm- 



Figure 1: Functions W q = W q ^ and Wcm = Wcm i n our bag model; e q = 
2.04-^, E = 8.17-^, in = 7.62-h, where R is the bag radius. 

Functions W q ± {x) and W^^(x) calculated using the bag model wavefunc¬ 
tions are shown in Fig. [!]. 


6 Consistency of the model—sum rules 

Taking Fourier transform of Eq. © we obtain 

T^ia) = J Tqm{ ax)W^ff (x) dx, (26) 
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where 


T^{a) = e“« a j d 3 rip q (r)^ 0 r^ q (r + e z a), (27) 

T CAi( a ) = & lEa J d 3 f^!CMCM{r + e z a). (28) 


For wavefunctions normalised to 1 we obtain T ( ^ LJ -(0) = 1 = 0), and hence 

from (|26| ) 

r t ±a -(o) 

j = 1 - < 29 > 

Since has an interpretation as the quark distribution function, the above 

equation stands for obvious normalisation condition (or from the experimental 
point of view Bjorken |b|] or Adler jl^] sum rule). It should be stressed that 
other methods of obtaining the structure functions from the bag model have 
some problems with fulfilling this conditon Q] . 

Differentiating (^g|) n times and putting a = 0 we obtain 

r jvi-xO-Dfo') 

/ W e ± f j(x)x n ~ 1 dx = q , (30) 

J Tqm 1 } ( 0 ) 


where T ±±< ' n ^ 
form 


' j_y 

k da J 


T ±a -. The above expression may be rewritten in the 


Mefj 


(n) = 



(31) 


where M ±A -(n) = J W ±± (x)x n ~ 1 dx (Mellin moments). It is well known that 
in the bag model without any improvements for the centre of mass motion one 
gets 2M^- ± (2) = This remains true in any independent particle model with 
the “volume” energy with energy-momentum tensor of the form: 


TV = b(r)g 

where b(r) is chosen so to guarantee conservation of the total (quarks + “vol¬ 
ume”) energy-momentum tensor 

d v (ZT£ v +T? , )= 0. (32) 

On the other hand energy-momentum tensor of our nucleon moving in the 
auxiliary potential cannot be conserved, and so corresponding 2 Mqm{2 ) ^ 1. 
But if we modify the condition ([l2|) by taking into account CM motion 

3,(3 T£ v + T# v )=d v T£ v M , (33) 
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then we obtain 


M , J - L 


e " (2) M^( 2 ) 


M - 1 - 1 

q 


( 2 ) 


(34) 


although in such case 2M g L± (2) ^ 

On the other hand, equation ( p3| ) together with (25) allows us to find im¬ 
proved energy E of our bound nucleon. For N identical quarks it satisfies the 
following condition 


E — Er | — 


Ne a E 0 - Ne a 


E N 

where Eq is the nucleon energy calculated using (|ii 
For the bag model with massless quarks we obtain 

e 0 = 2.04-h E 0 = 4e fl = 8.17, 


— —— R 
where R is the bag radius. 


E = ?. 62 i, 


(35) 


(36) 


7 Extracting effective quark distributions 

To calculate effective quark distributions, we may use equation (O) or (|26), but 
since and W t j — are slowly decreasing functions (see Fig.^lf), Eq. (2^) is 

much more useful for numerical calculations. This equation is in fact Fredholm 
equation of the first kind, generally ill-posed — its solutions may not exist or 
may be not unique. The standard method of solving such problems is to convert 
them into variational ones with additional regularization term [Q. Following |^| 
we introduce a smoothing functional 


KW'if] = l 



2 

da T A 



where A > 0 is a regularization parameter. 

Now our goal is to find minimum of the functional T\. Since for physical 
solution > 0, we impose this condition as a constraint for our variational 

problem. 

We have solved our problem using procedure PTIPR from ||. Results are 
plotted in Fig. |]. 

In Fig. H we have plotted l.h.s. of the equation (^), which can be written 
in the shorthand form as 


Tq = Tcm * Weff 


(38) 


and its r.h.s. corresponding to solutions for A = 10 5 and A = 10 6 plotted in 

Fig. |. 
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Figure 2: Effective quark distributions q{x ) in the bag model found using vari¬ 
ational method with different regularization parameter A = 10 —5 ,10 6 ,10~ 7 ; 
e q = 2.04-^, E = 7.62-^, to = E\ For comparison there is also shown “naive” 
quark distribution q{x) = 2W^~ ± (x) for e q = 2.04-^, E = 8.17-^, to = 7.62-^. 


8 QCD evolution 

The obtained quark distributions in the proton cannot be compared directly 
with data. According to the usual procedure we assume them to be twist two 
contribution to the structure functions at small scale Q 2 = Qq , and then we 
evolve them to a higher scale. Results of the leading order (LO) evolution for 
F -2 to Q 2 = 4GeV 2 are plotted in the Fig. [|. There also is shown for comparison 
parametrisation of GRV j2|. The chosen initial scales are Qq = 0.23GeV 2 (the 
initial scale of GRV in the LO), and Qq = 0.35GeV 2 , which we found as giving 
better predictions in the small x region. 

It may be easily seen that evolution of our quark distributions from Qq = 
0.23GeV 2 gives too steep F 2 in the small x region, but it is possible to find Qq 
giving much better agreement with the data. On the other hand, our F 2 differs 
drastically from the data for x > 1CP 3 . We attempted to cure this problem 
taking into account target mass effects. 

9 Reinterpretation of the results in terms of the 
Politzer’s variable £ 

Let us start with the remark, that since GRV parametrisation was obtained 
by evolving (according to QCD) initial distributions given at small Q 2 = /i 2 , 
we don’t have to perform the evolution ourselves—it is enough to compare our 
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Figure 3: Comparison of T q and Tqm * kF e // for two solutions for W e ff plotted 
in Fig. 2 


predictions with that of GRV at some small scale larger then fj 2 . 

Described in the previous section procedure of finding quark distributions 
assumes that their entire dependence (with fixed Bjorken x) on Q 2 comes from 
QCD. But due to simple kinetic target mass effects it cannot be true for smaller 
Q 2 . It was found many years ago that much better scaling may be obtained 
in terms of the variable f related to the Bjorken x by 


£ 


1 + 


2 x 



4m%x 2 ’ 


(39) 


where m p stands for target (proton) mass. 

The connection between usual F 2 {x, Q 2 ) and the scaling function Fjf (£, Q 2 ), 
which for large Q 2 coincides with F ^ , (in that limit £ = x) is given by |n| : 


F 2 (x,Q 2 ) = 


(1 + Ax 2 m 2 /Q 2 ) 3 / 2 
m. . 

+6 p 


F 2 S (^Q 2 ) 


Q 2 (1 + 4:X 2 m 2 /Q 2 ) 2 


+12 


TO’ 


Q 4 (l+4x 2 m2/Q2)5/2 


J t deF 2 s (e,Q 2 ) ( 40 ) 

J 1 dt' f dCFiiC^Q 2 ), 


where £ is the function of x , see Eq.(^). 

F 2 (x,Q 2 ) obtained from the above expression with .F!f(£) = where 

F 2 ff (x) = xW^fj (x) is the effective bag structure function obtained in the 
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Figure 4: Results of the LO QCD evolution of the bag structure functions from 
Qq to Q 2 = 4GeV 2 , compared with GRV parametrisation at Q 2 = 4GeV 2 . 
Parameters are: A qcd = 200MeV, n/ = 4. 


section 7, compared with GRV parametrisation for the same Q 2 , is shown in 
the Fig. ||. 


10 Conclusions 

We have proposed here a method of finding hadron structures functions in in¬ 
dependent particle models. The method leads to the structure functions with 
the correct support and exactly fulfils the normalisation condition for valence 
quark distributions. 

The structure functions found for the bag model in the static cavity ap¬ 
proximation differs drastically from the “naive” ones, calculated without any 
improvement for the centre of mass motion. They disagree with the data, but 
it is possible to obtain rough agreement in the small x region, evolving them 
using perturbative QCD. 

Much better agreement with parametrisation of GRV may be obtained, tak¬ 
ing into account that structure functions scale for small Q 2 in terms of the 
Politzer’s variable £. 

It should be stressed that our disagreement with the large Q 2 data arises 
mainly from quite different behaviour of our valence quark distributions at low 
Q 2 , which are finite for x —> 0, and the ones used by GRV, which are singular in 
that limit (although their xq(x) —> 0). In fact, any independent particle model 
gives finite q v {x —> 0). To obtain singular q(x) one has probably to take into 
account Regge trajectories for final states in the deeply inelastic scattering. 
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Figure 5: Valence quark part of F 2 (x) calculated using Eq.(43) for Q 2 = 
0.35GeV 2 and Q 2 = 4GeV 2 compared with the parametrisation of GRV. 
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